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Abstract 

The  Federal  Emergency  Management  Agency  (FEMA),  Region  III  office, 
has  initiated  a  study  to  update  the  coastal  storm  surge  elevations  within 
the  states  of  Virginia,  Maryland,  and  Delaware,  and  the  District  of 
Columbia.  The  area  includes  the  Atlantic  Ocean,  the  Chesapeake  Bay, 
including  its  tributaries,  and  the  Delaware  Bay.  This  effort  is  one  of  the 
most  extensive  coastal  storm  surge  analyses  to  date,  encompassing  coastal 
floodplains  in  three  states  and  including  the  largest  estuary  in  the  world. 
The  study  will  replace  outdated  coastal  storm  surge  Stillwater  elevations 
for  all  Flood  Insurance  Studies  in  the  study  area,  and  serve  as  the  basis  for 
new  coastal  hazard  analysis  and  —  ultimately  —  updated  Flood  Insurance 
Rate  Maps  (FIRMs).  Study  efforts  were  initiated  in  August  of  2008,  and 
are  expected  to  conclude  in  2012. 

The  storm  surge  study  considers  both  tropical  storms  and  extratropical 
cyclones  for  determination  of  return  period  storm  surge  elevations.  The 
historical  record  of  events  in  the  region  has  been  used  to  reconstruct  wind 
and  pressure  fields  for  significant  extratropical  cyclone  events  and  to 
develop  a  synthetic  suite  of  tropical  storm  tracks  and  associated  parameters 
to  drive  the  models.  This  will  support  the  statistical  analysis  of  storm  surge 
return  period  elevations.  This  report,  the  third  of  three  reports  comprising 
the  required  Submittal  1  documentation,  describes  development  of  the  wind 
and  pressure  fields  required  for  FEMA  Region  III  storm  surge  modeling. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional  purposes. 
Citation  of  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial  products. 
All  product  names  and  trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this  report  are  not  to 
be  construed  as  an  official  Department  of  the  Army  position  unless  so  designated  by  other  authorized  documents. 
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Preface 

This  study  was  conducted  for  the  Federal  Emergency  Management  Agency 
(FEMA)  under  Project  HSFE03-06-X-0023,  “NFIP  Coastal  Storm  Surge 
Model  for  Region  III”  and  Project  HSFE03-09-X-1108,  “Phase  II  Coastal 
Storm  Surge  Model  for  FEMA  Region  III.”  The  FEMA  technical  monitor 
was  Robin  Danforth. 

The  work  was  performed  by  the  Coastal  Processes  Branch  (HF-C)  of  the 
Flood  and  Storm  Protection  Division  (HF),  US  Army  Engineer  Research 
and  Development  Center  -  Coastal  and  Hydraulics  Laboratory  (ERDC- 
CHL).  At  the  time  of  publication,  Mark  Gravens  was  Chief,  Coastal 
Processes  Branch;  Dr.  Ty  Wamsley  was  Chief,  Flood  and  Storm  Protection 
Division;  Dr.  Jeffrey  L.  Hanson  was  the  Project  Manager;  William  R. 

Curtis  was  Technical  Director  for  Flood  and  Coastal  Storm  Research  and 
Development.  The  Acting  Director  of  CHL  was  Jose  E.  Sanchez. 

COL  Kevin  J.  Wilson  was  the  Commander  of  ERDC,  and  Dr.  Jeffery  P. 
Holland  was  the  Director. 
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1  Overview 

The  Federal  Emergency  Management  Agency  (FEMA)  is  responsible  for 
preparing  Federal  Insurance  Rate  Maps  (FIRMs)  that  delineate  flood 
hazard  zones  in  coastal  areas  of  the  United  States.  Under  Task  Orders 
HSFE03-06-X-0023  and  HSFE03-09-X-1108,  the  US  Army  Corps  of 
Engineers  (USACE)  and  project  partners  are  assisting  FEMA  in  the  develop¬ 
ment  and  application  of  a  state-of-the-art  storm  surge  risk  assessment 
capability  for  the  FEMA  Region  III  domain,  which  includes  the  Delaware 
Bay,  Chesapeake  Bay,  District  of  Columbia,  Delaware-Maryland- Virginia 
Eastern  Shore,  Virginia  Beach,  and  all  tidal  tributaries  and  waterways 
connected  to  these  systems.  The  goal  is  to  develop  and  apply  a  complete 
end-to-end  modeling  system,  with  all  required  forcing  inputs,  for  updating 
the  floodplain  levels  for  coastal  and  inland  watershed  communities.  Key 
components  of  this  work  include: 

1.  developing  a  high-resolution  Digital  Elevation  Map  (DEM)  for  Region  III, 
and  converting  this  to  an  unstructured  modeling  grid,  with  up  to  50  m 
horizontal  resolution,  for  use  with  the  production  system; 

2.  defining  the  Region  III  storm  hazard  in  terms  of  historical  extratropical 
cyclones  and  synthetic  hurricane  parameters; 

3.  preparing  an  end-to-end  modeling  system  for  assessment  of  Region  III 
coastal  storm  surge  hazards; 

4.  verifying  model  accuracy  on  a  variety  of  reconstructed  tropical  and 
extratropical  cyclone  events; 

5.  applying  the  modeling  system  to  compute  the  10-,  50-,  100-,  and  500-  year 
floodplain  levels;  and 

6.  developing  a  database  with  GIS  tools  to  facilitate  archiving,  distribution, 
and  analysis  of  the  various  storm  surge  data  products. 

Under  the  direction  of  FEMA  Region  III  Program  Manager,  Robin 
Danforth,  USACE  assembled  a  multi-organization  partnership  to  meet  the 
Region  III  objectives.  Work  on  this  project  has  made  extensive  use  of  the 
capabilities  and  technology  developed  for  the  North  Carolina  Floodplain 
Mapping  Program  (NCFMP).  The  availability  of  the  NCFMP  storm  surge 
modeling  system  (Blanton  et  al.  2008)  has  resulted  in  a  significant  cost 
savings  for  FEMA  Region  III.  Experts  in  the  fields  of  coastal  storm-surge, 
wind-driven  waves,  Geospatial  Information  Systems  (GIS),  and  high- 
performance  computational  systems  have  worked  together  in  this  effort. 
The  project  partners  and  their  primary  roles  are  listed  in  Table  1-1. 
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Table  1-1.  Study  team. 


Organization 

Contacts 

Primary  Role(s) 

US  Army  Corps  of  Engineers 

Field  Research  Facility  (USACE-FRF) 

Jeff  Hanson 

Mike  Forte 

Heidi  Wadman 

Project  Manager 

DEM  Construction 

Model  Validations 

Applied  Research  Associates/Intra  Risk  (ARA) 

Peter  Vickery 

Dhiraj  Wadhera 

Simulated  Hurricanes 

Simulated  Hurricanes 

ARCADIS 

Hugh  Roberts 

John  Atkinson 

Shan  Zou 

Modeling  Mesh 

Modeling  Mesh 

Modeling  Mesh 

Elizabeth  City  State  University 

Jinchun  Yuan 

Web/GIS 

Oceanweather 

Vince  Cardone 

Andrew  Cox 

Wind  Field 
Reconstructions 

Renaissance  Computing  Institute  (RENCI) 

Brian  Blanton 

Lisa  Stillwell 

Kevin  Gamiel 

Modeling  System 

DEM  Construction 

Data  base/Web/G  IS 

University  of  North  Carolina-  Chapel  Hill  (UNC- 
CH) 

Rick  Luettich 

Science  Consultant 

US  Army  Corps  of  Engineers 

District  Offices  (NAP,  NAO,  NAB) 

Jason  Miller 

Paul  Moye 

Jared  Scott 

Bathy/Topo  Data 
Inventory 

In  addition  to  the  study  team,  a  Technical  Oversight  Group  provided 
guidance  and  input  to  all  project  phases.  This  group  included  members 
from  the  following  organizations: 

1.  Chesapeake  Bay  Research  Consortium 

2.  Delaware  Flood  Mitigation  Program 

3.  Federal  Emergency  Management  Agency  (FEMA) 

4.  Dewberry,  Inc. 

5.  North  Carolina  Floodplain  Mapping  Program 

6.  USACE  Engineer  Research  and  Development  Center. 

Two  prior  Submission  1  stand-alone  reports  (Submission  1.1  by  Forte  et  al. 
2011  and  Submission  1.2  by  Blanton  et  al.  2011)  provided  details  on  the 
construction  of  the  bathymetric/topographic  digital  elevation  model  and 
preparation  of  the  modeling  system  including  development  of  the  Advanced 
CIRCulation  (ADCIRC)  mesh.  This  report  (Submission  1.3)  provides  the 
final  set  of  documentation  required  for  Report  No.  1  -  Scoping  and  Data 
Review.  The  contents  of  each  Submission  are  listed  in  Table  1-2.  Guidelines 
for  study  conduct  and  documentation  appear  in  FEMA  (2007). 
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Table  1-2.  Contents  of  the  Submission  1  report. 


Submission 

Title 

Contents 

1.1 

FEMA  Region  III  Coastal  Storm  Surge  Analysis: 
Study  Area  and  Digital  Elevation  Model  (DEM) 

Project  Overview 

Study  Area 

DEM  Development 

1.2 

FEMA  Region  III  Coastal  Storm  Surge  Analysis: 
Computational  System 

Modeling  System 

Mesh  Development 

1.3 

FEMA  Region  III  Coastal  Storm  Surge  Analysis: 
Storm  Forcing 

Hurricane  Parameters 
Extratropical  Cyclones 

The  following  chapters  describe  the  development  of  a  synthetic  suite  of 
tropical  storm  tracks,  their  associated  parameters,  and  the  reconstruction  of 
wind  and  pressure  fields  for  significant  extratropical  cyclone  events  in  the 
FEMA  Region  III  coastal  domain. 
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2  Hurricane  Parameters 

The  Joint  Probability  Method  (JPM)  for  simulating  hurricane  risk  has  been 
used  in  one  form  or  another  since  the  late  1960s.  The  original  JPM  applica¬ 
tion,  while  not  called  JPM,  was  developed  by  Russell  (1968)  for  predicting 
wave  loads  on  offshore  structures  in  the  Gulf  of  Mexico.  The  approach  used 
by  Russell  was  a  full  Monte  Carlo  simulation  where  hurricanes  were 
modeled  using  straight-line  segments  with  wind  and  wave  fields  computed 
using  hurricane  wind  and  wave  models.  This  methodology  was  first 
introduced  because  the  number  of  historical  events  (hurricanes)  at  any  one 
location  is  insufficient  to  enable  standard  statistical  techniques  (such  as 
extreme  value  analyses)  to  estimate  flood  risk,  wave  height  risk,  wind  speed 
risk,  etc.  For  coastal  risk  assessment,  the  introduction  of  long  duration 
synthetic  tracks  that  mimic  the  behavior  of  hurricanes  while  they  are  off¬ 
shore  (and  generating  a  wave  field)  was  introduced  by  Resio  et  al.  (2007). 
Modeling  the  full  storm  track  of  a  hurricane  for  the  purpose  of  developing  a 
wind  hazard  curve  was  first  introduced  by  Vickery  et  al.  (2000).  These 
simulation  methodologies  both  model  the  correlations  between  storm 
intensity  (central  pressure)  and  radius  to  maximum  winds  ( RMW ).  Vickery 
et  al.  (2009a)  modeled  a  relationship  between  RMW  and  the  Holland  B 
(Holland  1980)  parameter,  whereas  Resio  et  al.  (2007)  treated  the  Holland 
B  parameter  as  being  deterministic. 

The  JPM  approach  is  a  simulation  methodology  that  relies  on  the  develop¬ 
ment  of  statistical  distributions  of  key  hurricane  input  variables  (central 
pressure,  radius  of  maximum  winds,  translation  speed,  and  heading)  and  a 
sampling  from  these  distributions  to  develop  model  hurricanes.  Heading  is 
the  direction  of  motion  of  the  storm  defined  as  clockwise  from  north.  A 
storm  with  a  heading  of  90  deg  is  headed  east,  and  one  with  a  heading  of  - 
90  deg  is  heading  west.  Heading  is  defined  from  -180  deg  to  +180  deg, 
such  that  storms  with  a  positive  heading  have  an  eastward  component  and 
storms  with  a  negative  heading  have  a  westward  component.  The  simula¬ 
tion  results  in  a  family  of  modeled  storms  that  preserve  the  relationships 
between  the  various  input  model  components,  but  provides  a  means  to 
model  the  effects  and  probabilities  of  storms  that  have  not  yet  occurred. 
For  this  study,  a  method  known  as  Optimum  Sampling  (JPM-OS)  will  be 
used,  which  seeks  to  reduce  the  number  of  required  simulated  storms. 
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Figures  2-1  through  2-3  show  the  tracks  of  historic  hurricanes  affecting  the 
study  region.  Of  the  few  hurricanes  affecting  the  study  region,  none  were 
intense  hurricanes.  In  the  absence  of  sufficient  historical  data,  the  approach 
taken  here  is  to  use  the  synthetic  hurricane  model  described  in  Vickery  et  al. 
(2009a)  as  a  surrogate  to  develop  distributions  of  central  pressure 
difference,  translation  speed,  and  storm  heading  for  storms  affecting  the 
study  area.  Thus,  a  hybrid  approach  was  used  that  synthesizes  the 
efficiencies  of  the  JPM  with  a  database  of  storm  statistics  generated  for  a 
100,000-year  timeframe.  The  storm  statistics  were  obtained  using  the 
synthetic  storm  modeling  approach  described  in  Vickery  et  al.  (2009b). 

2.1  Overview  of  storm  sampling  and  statistical  analysis  approach 

The  occurrence  of  hurricanes  in  the  FEMA  Region  III  domain  is  rare  —  on 
average,  once  every  seven  years  —  and  the  occurrence  of  landfalling  hurri¬ 
canes  —  which,  in  general,  have  potential  to  produce  the  greatest  storm 
surges  —  is  rarer  still.  Those  hurricanes  that  have  impacted  the  region 
during  the  past  100  years  or  so,  for  which  data  are  available,  represent  a 
very  small  subset  of  the  hurricanes  that  are  possible.  Vickery  et  al.  (2000, 
2009b)  present  a  method  for  defining  what  the  possibilities  of  other 
hurricanes  are  in  terms  of  tracks,  intensity,  size,  translation  speed,  and 
radial  changes  in  the  wind  field,  based  on  the  historical  data  record  for  the 
larger  Atlantic  coast  region.  This  method  involved  stochastic  simulation  of 
100,000  hypothetical  hurricanes  using  a  Monte  Carlo  technique;  it  repli¬ 
cated  the  statistics  of  Atlantic  coast  hurricanes  that  have  influenced  the 
region  in  terms  of  tracks  and  other  hurricane  parameters. 

The  approach  that  is  being  adopted  in  this  study  for  characterizing  storm 
water  levels  and  waves  utilizes  very  high  resolution  computer  modeling  to 
maximize  accuracy  and  spatial  granularity  of  model  output.  Simulation  of 
100,000  storms  using  high  resolution  modeling  is  not  computationally 
feasible  and  would  require  several  years  to  complete,  so  the  number  of 
storms  had  to  be  reduced.  The  JPM  was  adopted  to  reduce  the  number  of 
storms.  The  JPM  is  a  simulation  methodology  that  relies  on  the 
development  of  statistical  distributions  of  key  hurricane  input  parameters 
(including  central  pressure,  radius  of  maximum  winds,  translation  speed 
and  headings  or  paths);  the  definition  of  discrete  values  for  these 
parameters  to  adequately  represent  their  statistical  distributions;  and  the 
definition  of  a  set  of  storms  that  have  a  unique  set  of  parameter  values  and 
that  are  specified  to  travel  along  a  discrete  idealized  path.  A  reduction  in  the 
number  of  storms  must  be  carefully  calculated  to  preserve  the  inherent 
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Figure  2-1.  Tracks  of  all  hurricanes  affecting  FEMA  Region  III  for  the  1940- 

2007  period. 


Figure  2-2.  Tracks  of  Category  2  hurricanes  affecting  FEMA  Region  III  for 
the  1940-2007  period. 
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Figure  2-3.  Tracks  of  Category  3  and  higher  hurricanes  affecting  FEMA 
Region  III  for  the  1940-2007  period. 


probabilities  of  hurricane  characteristics.  The  statistics  for  hurricane 
parameters  and  tracks  were  computed  using  the  method  of  Vickery  et  al. 
(2009b).  These  statistics  were  used  as  the  preservation  target  in  the  storm 
reduction  process.  Another  method  used  to  help  guide  the  number 
reduction  process  was  a  coarse  grid  storm  surge  model;  this  model 
simulated  50,000  of  the  100,000  storms.  Those  results  were  then  used  to 
estimate  exceedance  probabilities  for  storm  surges  throughout  the  Region 
III  domain.  The  team  eventually  reduced  the  JPM  set  to  156  storms  to 
preserve  the  statistics  of  both  the  hurricane  parameters  and  their  paths,  and 
the  storm  surge.  This  storm  set  reduction  process  is  called  an  optimal 
sampling  (OS)  approach,  or  JPM-OS. 

Through  an  iterative  process,  the  50,000  storms  were  reduced  to  17,644; 
that  number  was  further  reduced  to  a  JPM-OS  set  of  468  storms  —  each 
storm  with  a  unique  pathway  and  set  of  hurricane  input  parameters.  The 
468  storms  consisted  of  three  classes  of  storms:  one  class  that  made  landfall 
in  North  Carolina;  one  class  that  made  landfall  along  the  Virginia/ 
Maryland/ Delaware  coast;  and  a  third  class  —  bypassing  hurricanes  —  that 
passed  just  offshore  and  near  the  coast  but  did  not  make  landfall.  The 
468-storm  set  included  storms  that  made  landfall  either  at  evenly  spaced 
intervals  along  the  two  coastline  segments,  or  —  in  the  case  of  bypassing 
hurricanes— at  equally  spaced  points  along  a  shore-perpendicular  transect. 
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An  additional  reduction  in  the  number  of  storms  (to  156)  was  made  by 
eliminating  variability  in  the  hurricane  parameter  (the  Holland  B 
parameter)  that  controls  radial  wind  speed  changes.  Eliminating  this 
variability  had  little  effect  on  storm  surge  probabilities.  This  reduction  left 
the  storms  at  unevenly  spaced  intervals  for  both  landfalling  and  bypassing 
hurricanes.  Analysis  was  done  to  examine  the  effect  of  using  unevenly 
spaced  storms,  compared  to  evenly  spaced  storms,  and  the  differences  in 
terms  of  storm  surge  probabilities  were  found  to  be  negligible.  The  final 
JPM-OS  storm  set  consisted  of  156  unequally  spaced  storms.  The  decision 
to  adopt  unevenly  spaced  storms  was  also  done  to  avoid  rerunning  storms 
through  the  high  resolution  storm  surge  and  wave  models. 

The  details  of  the  progression  to  the  final  JPM-OS  set  of  156  storms  — 
including  the  results  of  analyses  performed  to  inform  the  storm  reduction 
process  —  are  described  in  more  detail  in  the  following  sections. 

2.2  Stochastic  tracks 

The  approach  for  conducting  the  hurricane  simulations  for  storms  affecting 
the  Chesapeake  Bay  area  and  coastal  Virginia  and  Delaware  makes  use  of 
the  simulated  set  of  hurricane  tracks  as  given  in  Vickery  et  al.  (2009b).  The 
simulation  methodology  is  outlined  in  Figure  2-4.  Using  this  methodology, 
storms  are  initiated  in  the  Atlantic  Basin  at  the  locations  and  times  observed 
in  the  historical  record.  The  storms  are  moved  across  the  ocean  using  a 
Markov-type  model,  where  the  position,  speed,  and  heading  of  each  storm 
at  the  next  time  step  is  determined  by  sampling  from  distributions  of 
changes  in  speed  and  heading  that  are  primarily  functions  of  the  speed  and 
heading  of  each  storm  at  its  respective,  current  location.  Different  distribu¬ 
tions  of  changes  in  storm  speed  and  heading  are  used  in  different  five 
degree  squares  encompassing  the  Atlantic  Basin.  The  storm  intensity,  as 
defined  using  central  pressure,  is  limited  by  potential  intensity  theory 
(Emanuel  1988). 

A  simple  ocean  mixing  model  (Emanuel  et  al.  2004)  is  used  to  mix  cooler 
waters  from  beneath  the  surface  of  the  ocean  to  limit  intensity  changes. 

Along  the  coastline  of  the  United  States,  the  model  was  validated  through 
comparisons  of  statistical  distributions  of  storm  heading,  translation  speed, 
distance  parameters,  central  pressure,  and  frequency.  These  validations 
were  performed  using  both  landfalling  hurricane  data  as  well  as  the  storm 
data  computed  at  the  time  a  storm  is  closest  to  the  center  of  a  circle  having 
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Figure  2-4.  Overview  of  synthetic  storm  simulation  methodology  described  in 
Vickery  et  al.  (2009b).  In  this  figure,  /  represents  relative  intensity  as  defined  in 
Emanuel  (1988),  RMW  is  the  radius  to  maximum  winds,  Pc  is  the  central  pressure, 

B  is  the  Holland  (1980)  pressure  profile  parameter,  0  is  the  storm  heading,  and  Vt 
is  the  translation  speed  of  the  storm. 

a  radius  of  250  km.  These  sample  circles  were  located  at  approximately 
50  nautical  mile  increments  along  the  entire  coastline  of  the  United  States. 
This  synthetic  storm  model  was  shown  to  reproduce  the  statistics  of 
hurricanes  affecting  the  U.S.  coastline  from  1900  through  2008,  and  is 
capable  of  simulating  hundreds  of  thousands  of  years  of  storms  whose 
statistical  characteristics  match  those  of  the  historic  record. 

These  simulated  hurricanes,  whose  statistical  properties  are  consistent  with 
those  affecting  the  region,  were  used  instead  of  the  extremely  limited  data 
associated  with  historical  hurricanes  affecting  the  study  area.  Figure  2-5 
shows  milepost  locations  (MP  150  to  MP  3100)  and  26  sections  of  the  U.S. 
coastline.  This  study  was  influenced  by  the  statistics  of  hurricanes  affecting 
MPs  2100  through  2400  and  coastline  segments  21,  22,  23,  and  24.  Figure 
2-6  shows  the  areas  covered  by  the  250  km  radius  circles  centered  on  the 
mileposts  in  the  region. 
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Figure  2-5.  Milepost  and  coastal  segments.  Each  coastal  segment  (lower  map)  is 
about  200  km  in  length  and  contains  between  one  and  two  milepost  markers, 
depending  on  the  exact  location  of  the  segment. 


Figure  2-7  presents  a  comparison  of  modeled  (Vickery  et  al.  2009b)  and 
observed  central  pressures  from  hurricanes  making  landfall  along  each  of 
the  26  coastal  segments  shown  in  Figure  2-5. 

As  indicated  in  Figure  2-7,  only  four  hurricanes  (with  central  pressures 
less  than  990  mbar)  have  made  landfall  along  coastline  Segment  22,  and 
no  hurricanes  made  landfall  along  coastline  Segment  23.  The  hurricane 
simulation  model,  which  has  been  run  for  a  period  of  100,000  years, 
generates  synthetic  hurricanes  that  make  landfall  along  these  coastal 
segments.  The  relatively  large  number  of  simulated  years  (100,000)  was 
originally  required  to  enable  estimates  of  long  return  period  (-2000  years 
or  longer)  wind  speeds. 
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Figure  2-6.  Close-up  of  areas  encompassed  by  the  250  km  radius  sub- 
regions  used  to  characterize  storms  affecting  the  study  region.  Darker 
coastline  defines  the  FEMA  Region  III  coastline. 


Figure  2-8  presents  a  comparison  of  modeled  and  observed  pressures  at  the 
time  the  storms  are  closest  to  the  indicated  milepost.  The  comparisons  also 
show  good  agreement  between  the  modeled  and  historical  central  pressure 
differences.  Unlike  the  comparisons  shown  in  Figure  2-7,  which  include 
landfalling  hurricanes  only,  the  model  and  historical  data  shown  in 
Figure  2-8  includes  all  tropical  cyclones.  Hence,  the  model  provides  a  much 
larger  data  set  of  historic  storms  to  use  for  comparisons  of  modeled  and 
historical  pressure  data.  Tests  for  equivalence  of  means  and  variance  of  the 
synthetic  and  historical  pressures  shown  in  Figure  2-8  indicate  that  there  is 
no  statistically  significant  difference  (evaluated  at  the  95th  percent 
confidence  interval)  between  the  two  data  sets.  The  two  distributions  at  all 
mileposts  pass  the  Kolmogorov-Smirnov  test  for  equivalence  at  the  95th 
percent  confidence  level. 
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Figure  2-7.  Comparison  of  modeled  and  observed  (1900-2006)  central  pressures  at 
landfall  vs.  return  period  along  ~200km  coastal  segments.  Gray  lines  show  90  % 
confidence  range  derived  from  the  modeled  empirical  distribution.  If  the  upper  gray  line 
is  not  shown,  the  five  percent  (lower  bound)  simulated  pressures  are  less  than  990  hPa. 

Figures  2-9  and  2-10  present  comparisons  of  the  statistics  (translation 
speed  and  heading)  of  all  simulated  and  observed  tropical  cyclones  passing 
within  250  km  of  Mileposts  2250  and  2300  near  the  North  Carolina  and 
Virginia  border.  As  indicated  in  Figures  2-7  through  2-9,  the  statistical 
distributions  of  the  key  tropical  cyclone  parameters  derived  from  the 
100,000  year  simulation  are  seen  to  be  statistically  equivalent  to  those  of 
the  short  106  year  historical  record,  but  because  the  synthetic  storm  set  was 
developed  using  a  100,000  year  simulation,  there  is  a  relatively  large 
number  of  landfalling  and  close  by-passing  hurricanes  that  can  be  used  to 
define  the  statistical  distributions  of  the  key  hurricane  parameters  for  use  in 
the  JPM  models. 
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Figure  2-8.  Cumulative  distributions  of  central  pressures  from  simulated  and  historical 
tropical  cyclones  passing  within  250  km  of  the  indicated  milepost.  Central  pressure 
represents  the  value  at  the  time  a  storm  is  closest  to  the  indicated  milepost. 
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Figure  2-9.  Comparisons  of  translation  speeds  from  simulated  and  historical  tropical  cyclones 
passing  within  250  km  of  the  indicated  milepost.  Dashed  lines  show  the  bounds  associated 

with  the  Kolmogorov-Smirnov  test. 
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Figure  2-10.  Comparisons  of  storm  heading  from  simulated  and  historical  tropical  cyclones 
passing  within  250  km  of  the  indicated  milepost.  Dashed  lines  show  the  bounds  associated 
with  the  Kolmogorov-Smirnov  test.  Heading  is  defined  as  the  direction  of  storm  motion 
measured  clockwise  from  North.  Storms  with  a  positive  heading  have  an  eastward  component. 


In  this  study,  the  simulated  storms  are  used  to  define  the  statistical 
distributions  and  frequencies  of  hurricanes  affecting  the  Chesapeake  Bay 
and  coastal  Virginia  and  Delaware  regions.  Using  the  synthetic  storm 
model,  100,000  years  of  storm  tracks  were  generated  and  provided  to  the 
Renaissance  Computing  Institute  (RENCI)  to  be  used  in  conjunction  with 
a  coarse  grid  in  the  Advanced  CIRCulation  storm  surge  model  (ADCIRC). 
The  storm  track  data  was  also  used  in  a  simple  Holland  (1980)  windfield 
model  to  generate  hurricane-induced  water  elevations  in  the  study  region. 

2.3  Coarse  grid  results 

A  coarse  ADCIRC  grid  was  used  to  provide  a  baseline  from  which  to 
compare  JPM  subset  simulations.  The  coarse  grid  was  derived  from  the 
standard  ADCIRC  east  coast  2001  grid  by  extracting  the  nodes  and  elements 
in  the  Region  III  and  adjacent  areas.  There  are  8178  nodes  and  15450 
elements,  with  a  minimum  depth  of  0.5  m.  The  grid  extends  northward  to 
the  Delaware/New  Jersey  coast,  southward  to  Cape  Hatteras,  and  offshore 
past  the  continental  shelf  break,  and  includes  only  open  water. 
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The  original  storm  set  provided  to  RENCI  for  use  with  the  coarse  ADCIRC 
grid  contained  over  57,000  storms;  17,644  produced  a  storm  surge  of  at 
least  one  meter  somewhere  in  the  region  “P”  shown  in  Figure  2-11. 

Figure  2-12  presents  the  storm  surge  elevation  as  a  function  of  the  return 
period  at  each  of  the  five  points  shown  in  Figure  2-11,  as  well  as  the  maxi¬ 
mum  storm  surge  within  the  region  “P,”  computed  from  the  full  storm  set 
(-57,000  storms)  in  the  upper  plot  and  the  reduced  storm  set 
(17,644  storms)  in  the  lower  plot.  The  full  storm  set  includes  all  simulated 
tropical  cyclones  that  pass  within  250  km  of  any  of  the  mileposts  shown  in 
Figure  2-6.  A  comparison  of  the  storm  surge  hazard  curves  resulting  from 
the  full  and  reduced  storm  sets  indicates  that  the  elimination  of  nearly 
40,000  weaker  or  distant  storms  has  virtually  no  effect  on  the  storm  surge 
elevations  for  return  periods  of  10  years  and  longer. 

The  storms  affecting  the  region  were  divided  into  three  classes:  namely, 
storms  making  landfall  along  the  North  Carolina  coast,  storms  making  land¬ 
fall  along  the  Virginia/Delaware  coast,  and  by-passing  storms.  The  line  seg¬ 
ments  associated  with  these  three  classes  of  storms  are  shown  in  Figure  2-13. 


Figure  2-11.  Region  used  for  the  coarse  grid  ADCIRC  run. 
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Figure  2-12.  Storm  surge  elevation  vs.  return  period  from  the  coarse  grid  ADCIRC  run. 
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Figure  2-13.  Coastal  segments  used  to  define  hurricane  tracks  for  the 
Region  III  Flood  Study. 


2.4  JPM  track  selection 

Hurricane  track  parameters  (central  pressure,  translation  speed,  and 
heading)  and  their  weights  were  developed  for  the  three  line  segments  used 
to  define  the  hurricane  risk  for  the  study.  Along  each  of  the  two  coastal 
segments,  the  statistical  distributions  of  heading  and  central  pressure  are 
approximated  using  three  discrete  values.  The  distribution  for  translation 
speed  is  characterized  using  two  discrete  values.  In  the  first  sets  of  simula¬ 
tions  performed  during  the  2010  timeframe,  the  distributions  of  the  radius 
to  maximum  winds  and  the  Holland  B  parameter  were  defined  using  the 
mean  value  and  the  mean  ±1.2205,  where  ob  is  the  error  term  from  the 
regression  equations  relating  B  to  RMW  and  latitude  is  as  described  in 
Vickery  and  Wadhera  (2008).  Using  the  mean  ±1.2205  yields  a  set  of  B 
values  that  matches  the  mean  and  the  variance  of  the  error  term  from  the 
regression  model. 

The  Holland  B  parameter  is  modeled  using  Equation  24  in  Vickery  and 
Wadhera  (2008)  and  the  radius  to  maximum  winds  is  modeled  using 
Equation  11  in  Vickery  and  Wadhera  (2008).  Figure  2-14  presents  the 
mean  and  ±1.220  bounds  for  the  Holland  B  model  plotted  vs.  storm  radius 
to  maximum  winds.  Also  provided  in  Figure  2-14  are  the  historical  values 
of  B  obtained  from  hurricane  model  hindcasts,  as  described  in  Vickery  and 
Wadhera  (2008). 
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RMW  (km)  | 

Figure  2-14.  Modeled  and  observed  values  of  the  Holland  B  parameter  plotted  vs.  radius  to 
maximum  winds.  Historical  data  from  recent  Atlantic  coast  landfalling  and  by-passing 

hurricanes  (excludes  Florida). 

These  B  values  were  used  in  the  hurricane  wind  field  model  described  in 
Vickery  et  al.  (2009a),  which  is  the  same  model  used  in  the  FEMA  Region 
III  production  runs,  as  indicated  in  the  FEMA  Region  III  report  entitled 
Coastal  Storm  Surge  Analysis:  Computational  System,  Report  2: 
Intermediate  Submission  1.2,  February  2011.  As  indicated  in  Figure  2-14, 
both  the  model  and  the  data  indicate  that  B  decreases  with  increasing 
radius.  It  is  also  noteworthy  that  the  limits  defined  by  ±1.22 ob  match  the 
limits  of  the  limited  historical  data. 

Figure  2-15  presents  the  mean  and  ±1.220  bounds  for  the  RMW  model 
plotted  vs.  central  pressure.  Also  provided  in  Figure  2-15  are  the  historical 
values  of  RMW  obtained  from  hurricane  model  hindcasts,  as  described  in 
Vickery  and  Wadhera  (2008).  Both  the  model  and  the  data  indicate  that 
RMW  decreases  with  decreasing  central  pressure.  RMW  also  increases 
with  increasing  latitude.  As  in  the  case  of  B,  it  was  noted  that  the  limits 
defined  by  ±1.220  match  the  limits  of  the  limited  historical  data. 

This  combination  of  storm  parameters  (three  radii  and  three  B  models) 
results  in  a  total  of  nine  simulated  hurricanes  for  each  pressure-heading- 
translation  speed  combination.  In  all  cases,  the  weights  were  first  chosen  so 
that  the  means  and  standard  deviations  of  the  weighted  discrete  distribu¬ 
tions  matched  those  of  the  full  continuous  distributions  derived  from  the 
stochastic  simulation  model.  An  acceptable  set  of  468  synthetic  storms  was 
developed  using  this  approach. 
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Figure  2-15.  Modeled  and  observed  values  of  radius  to  maximum  plotted  vs.  central 
pressures.  Historical  data  from  recent  Atlantic  coast  landfalling  and  by-passing  hurricanes 

(excludes  Florida). 


Near  the  end  of  the  track  development  stage,  a  sensitivity  study  was 
performed  where  2/3  of  tracks  were  eliminated  (those  corresponding  to 
±1.2205)  —  leaving  only  the  mean  regression  model  for  B  —  relating  B  to 
RMW and  latitude.  The  resulting  storm  surge  estimates  associated  with  a 
return  period  of  100  years  differed  little  from  those  developed  using  the 
full  set,  and  a  decision  was  made  to  proceed  with  this  reduced  set  of 
synthetic  hurricane  tracks.  This  reduced  set  of  156  storms  contains  1/3  the 
number  of  the  full  468  storm  JPM  set. 

As  will  be  discussed  in  Section  2.5,  the  storm  surge  estimates  resulting 
from  the  reduced  storm  set  used  in  the  JPM  modeling  (156  individual 
hurricanes)  are  verified  through  comparisons  to  storm  surge  results 
derived  from  17,644  simulated  storms  (this  number  was  a  subset  of  the  full 
Monte  Carlo  simulation  storm  set). 

Figures  2-16  and  2-17  present  plots  of  the  maximum  storm  surges  produced 
using  the  coarse  ADCIRC  model  runs  plotted  as  a  function  of  each  of  the 
three  primary  variables  (heading,  translation  speed,  and  central  pressure 
difference),  for  each  of  the  two  coastal  landfalling  segments.  Also  shown  in 
Figures  2-16  and  2-17  are  the  discrete  and  cumulative  distribution  functions 
describing  storm  heading,  central  pressure  difference,  and  translation 
speed.  The  JPM  representation  of  the  Virginia/Delaware/New  Jersey 
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Figure  2-16.  Plots  of  storm  surge  vs.  hurricane  parameters  and  comparisons  of  continuous 
and  discrete  distributions  of  storm  heading,  central  pressure  difference  (c/P),  and  translation 
speed  (c)  for  Virginia,  Delaware  and  New  Jersey  landfalling  hurricanes.  Heading  is  defined  as 
the  direction  of  motion  of  the  storm  measured  clockwise  from  North.  Storms  with  a  positive 
heading  have  an  eastward  component.  A  total  of  four  JPM  representations  of  the  translation 
speed  are  used  to  model  the  translation  speed  for  the  five  different  headings.  The  same 
translation  speed  distribution  is  used  for  headings  -45  and  -30  deg. 
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Figure  2-17.  Plots  of  storm  surge  vs.  hurricane  parameters  and  comparisons  of  continuous 
and  discrete  distributions  of  storm  heading,  central  pressure  difference  and  translation  speed 
for  North  Carolina  landfalling  hurricanes.  Heading  is  defined  as  the  direction  of  motion  of  the 
storm  measured  clockwise  from  North.  Storms  with  a  positive  heading  have  an  eastward 
component.  A  total  of  three  JPM  representations  of  the  translation  speed  are  used  to  model 
the  translation  speed  for  the  three  different  headings. 


landfalling  storms  is  seen  to  have  a  distinct  bias  towards  more  intense 
(greater  central  pressure  differences  or  lower  central  pressures)  landfalling 
storms  than  the  results  of  the  stochastic  set.  This  distinct  bias  is  intentional 
and  was  required  to  correct  for  a  negative  bias  in  the  resulting  JPM  storm 
surge  estimates  (based  on  the  full  storm  set  runs  described  in  Chapter  2.2). 
In  the  North  Carolina  case,  no  central  pressure  adjustment  was  required. 
Figure  2-17  shows  that  the  JPM  tracks  reproduce  a  weak  correlation 
between  heading  and  central  pressure  that  was  not  evident  along  other 
coastal  segments.  Heading  is  defined  as  the  direction  of  storm  motion 
measured  clockwise  from  North,  and  storms  with  an  eastward  component 
in  their  motion  have  positive  heading  values. 
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Figure  2-16  (the  Virginia/Delaware/New  Jersey  landfalling  storms) 
demonstrates  that  there  is  a  different  translation  speed  probability 
distribution  used  for  each  heading  (the  translation  speed  distributions  for 
headings  -30  and  -45  deg  are  the  same),  resulting  in  a  set  of  tracks  that 
maintains  the  heading-translation  speed  correlation  evident  in  the 
stochastic  modeling  and  in  the  historical  record.  The  r2  associated  with  a 
linear  model  describing  the  relationship  between  heading  and  translation 
speed  is  0.16;  consequently,  the  team  decided  that  this  correlation  should  be 
retained  in  the  JPM  representation  of  the  hurricanes.  In  the  case  of  the 
North  Carolina  landfalling  hurricanes,  the  corresponding  r2  is  only  0.06, 
and  any  relationship  between  heading  and  translation  speed  was  ignored. 

To  reproduce  both  the  mean  and  the  variance  of  the  translation  speed  for 
the  Virginia/ Delaware/New  Jersey  storms  using  only  two  discrete  values, 
the  separation  between  the  lowest  and  highest  translation  speed  was  found 
to  be  relatively  large  for  the  eastward  heading  storms.  This  large  separa¬ 
tion  is  required,  as  both  the  mean  and  the  variance  of  the  translation 
speeds  for  these  eastward  heading  (re-curving)  storms  are  higher  than 
those  for  storms  having  a  more  westward  heading. 

In  the  case  of  bypassing  hurricanes,  the  distributions  of  the  hurricane’s 
parameters  are  characterized  using  only  one  heading,  two  values  of  central 
pressure,  and  one  translation  speed. 

Figure  2-18  presents  plots  of  the  maximum  storm  surges  produced  using 
the  simple  ADCIRC  model  runs  performed  earlier  plotted  as  a  function  of 
each  of  the  three  primary  variables  (heading,  translation  speed,  and 
central  pressure),  for  the  by-passing  segment.  This  figure  also  shows  the 
discrete  and  continuous  cumulative  distribution  functions  describing 
storm  heading,  central  pressure  difference,  and  translation  speed. 

Figure  2-19  shows  the  cumulative  distribution  functions  (CDFs)  of  the 
landfall  location  (or  coast  crossing  position)  of  the  hurricane  tracks  along 
each  of  the  three  coastal  segments.  Three  different  storm  position  models 
are  used  to  describe  the  track  landfall  location  along  the  North  Carolina 
and  Virginia  coastal  segments,  but  only  one  model  is  used  to  define  the 
positioning  of  the  bypassing  hurricanes.  The  track  spacing  in  all  cases  is 
obtained  using  equally  spaced  positions  along  the  ordinate  in  Figure  2-19, 
which  yields  non-equally  spaced  positions  along  the  coastal  line  segment. 
If  the  underlying  distributions  of  the  track  coast  crossings  were  uniform, 
this  equal  spacing  along  the  ordinate  approach  would  lead  to  uniform 
spacing  along  the  length  of  the  coastal  segments. 
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Figure  2-18.  Plots  of  storm  surge  vs.  hurricane  parameters  and  comparisons  of  continuous  and 
discrete  distributions  of  storm  heading,  central  pressure  difference  and  translation  speed  for  by¬ 
passing  hurricanes.  Heading  is  defined  as  the  direction  of  motion  of  the  storm  measured 
clockwise  from  North.  Storms  with  a  positive  heading  have  an  eastward  component. 


The  numeric  values  of  the  modeled  storm  parameters  and  their  statistical 
weights  are  given  in  Tables  2-1  and  2-2  for  the  landfalling  coastal  segments. 
The  parameter  values  and  weights  for  the  by-passing  hurricanes  are  given  in 
Table  2-3. 

The  along  track  variation  of  central  pressure  in  the  JPM  hurricanes  was 
treated  using  the  pre-landfall  central  pressure  difference  variation  model 
developed  during  the  North  Carolina  flood  study.  Once  a  hurricane  makes 
landfall,  the  central  pressure  difference  decreases  as  the  hurricane  moves 
inland.  The  post-landfall  filling  (increase  in  central  pressure  and  RMW 
and  a  decrease  in  B)  was  modeled  with  the  filling  models  given  in  Vickery 
(2005).  In  the  implementation  of  the  filling  model  for  this  study,  the 
random  error  term  in  the  Vickery  (2005)  model  was  included  rather  than 
simply  using  the  mean  filling  rate  model. 
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Figure  2-19.  Cumulative  distribution  functions  of  tracks  landfalling  along  the  three 
coastal  segments.  The  zero  position  is  the  leftmost  portion  of  each  line  segment  shown 

in  Figure  2-13. 


Table  2-1.  Hurricane  parameter  values  and  weights  for  Virginia,  Delaware,  and  New  Jersey 

landfalling  hurricanes. 


Heading  Degrees 
CW  from  North 

Weight 

Central  Pressure 
Difference  (mbar) 

Weight 

Trans  Speed 
(m/sec) 

Weight 

-75 

0.09 

34,  51,  65 

0.44,  0.35,  0.21 

4.0,  7.2 

0.60,  0.40 

-45 

0.06 

34,  51,  65 

0.44,  0.35,  0.21 

4.7,  9.5 

0.65,  0.35 

-30 

0.12 

34,  51,  65 

0.44,  0.35,  0.21 

4.7,  9.5 

0.65,  0.35 

-10 

0.27 

34,  51,  65 

0.44,  0.35,  0.21 

4.8,  10.5 

0.57,  0.43 

-15 

0.47 

34,  51,  65 

0.44,  0.35,  0.21 

5.0,  13.0 

0.45,  0.55 

Table  2-2.  Hurricane  parameter  values  and  weights  for  North  Carolina  landfalling  hurricanes. 


Heading  Degrees 
CW  from  North 

Weight 

Central  Pressure 
Difference  (mbar) 

Weight 

Trans  Speed 
(m/sec) 

Weight 

-35 

0.37 

38,  56,  75 

0.52,  0.40,  0.08 

4.6,  11.5 

0.66,  0.34 

0 

0.20 

38,  56,  75 

0.47,  0.42,  0.11 

4.6,  11.5 

0.66,  0.34 

22 

0.43 

38,  56,  75 

0.38,  0.50,  0.12 

4.6,  11.5 

0.66,  0.34 

Table  2-3.  Hurricane  parameter  values  and  weights  for  by-passing  hurricanes. 


Heading  Degrees 
CW  from  North 

Weight 

Central  Pressure 
Difference  (mbar) 

Weight 

Trans  Speed 
(m/sec) 

Weight 

12 

1 

42,67 

0.75,0.25 

4.0,11.4 

0.65,0.35 
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Figures  2-20  and  2-21  present  comparisons  of  probability  distributions  of 
heading,  central  pressure,  etc.,  at  the  time  of  landfall  along  the  North 
Carolina  and  Virginia/Delaware/New  Jersey  coastal  segments.  This 
information  is  the  result  of  both  the  stochastic  and  the  JPM  models.  The 
comparisons  show  that  the  limited  numbers  of  JPM  storms  combined  with 
the  weights  reproduce  the  CDFs  from  the  stochastic  simulation  reasonably 
well,  with  the  exception  that  the  full  variance  of  B  is  not  reproduced,  since 
B  is  modeled  using  the  mean  regression  model  only.  Comparisons 
resulting  from  the  468  storm  set  are  provided  in  Appendix  C,  where  the 
full  distribution  of  B  is  well  reproduced,  with  the  addition  of  the  other  316 
storms  required  when  modeling  the  mean  ±1.2205. 


-180  -135  -90  -45  0  45  90  135  180 

Storm  Heading 


B 


0 


20  40dp(mb)60 


80  100 


Figure  2-20.  Comparison  of  JPM  and  stochastic  representations  of  key  hurricane  parameters 
at  landfall  along  the  VA/DE/NJ  coastal  segment.  Note  that  although  only  the  mean 
representation  of  the  Holland  B  parameter  was  modeled  in  the  JPM  set,  some  variation  in  B  is 
observed  due  to  the  correlation  of  B  with  RMW  and  latitude.  Heading  is  defined  as  the 
direction  of  storm  motion  measured  clockwise  from  North.  Storms  with  a  positive  heading 
have  an  eastward  component.  The  translation  speed  is  denoted  c,  the  central  pressure 
difference  is  denoted  c/P,  and  RMW  is  the  radius  to  maximum  winds. 
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Figure  2-21.  Comparison  of  JPM  and  stochastic  representations  of  hurricane  parameters  at 
landfall  along  the  NC  coastal  segment.  Note  that  although  only  the  mean  representation  of 
the  Holland  B  parameter  was  modeled  in  the  JPM  set,  some  variation  in  B  is  observed  due  to 
the  correlation  of  B  with  RMW  and  latitude.  Heading  is  defined  as  the  direction  of  storm 
motion  measured  clockwise  from  North.  Storms  with  a  positive  heading  have  an  eastward 
component.  The  translation  speed  is  denoted  c,  the  central  pressure  difference  is  denoted 
dP,  and  RMW  is  the  radius  to  maximum  winds. 


Recall  that  the  JPM  central  pressures  shown  in  Figure  2-20  are  shifted 
toward  being  more  intense  than  the  stochastic  central  pressures,  which 
were  found  through  a  trial  and  error  process  to  be  necessary  to  correct  a 
negative  bias  in  the  100-year  return  period  results. 

2.5  Comparison  of  storm  surge  hazard  curves  derived  using  the  JPM 
approach  to  those  derived  from  a  full  stochastic  storm  set 

Results  presented  in  this  section  are  for  a  set  of  156  storm  tracks,  which 
were  originally  a  subset  of  the  base  (or  full  JPM)  storm  set  containing  468 
simulated  tracks.  In  the  full  JPM  storm  set,  the  variation  of  the  Holland  B 
parameters  was  represented  by  modeling  the  mean  value  and  the  mean 
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±1.22  standard  deviations  but  —  here  in  this  reduced  set  —  only  the  mean 
representation  of  the  Holland  B  parameter  is  maintained.  In  the  selection 
of  the  156  storms,  the  tracks  associated  with  the  upper  and  lower  bound 
representations  of  B  were  removed. 

Figure  2-22  presents  the  nine  sets  of  tracks  used  in  this  156  storm  set.  To 
verify  that  the  weighted  simulation  methodology  yields  reasonable  results, 
and  to  determine  the  return  period  range  over  which  the  results  can  be 
considered  valid,  the  team  compared  the  estimates  of  the  predicted  hurri¬ 
cane-induced  storm  surge  (derived  using  the  JPM  model  hurricanes)  to 
those  predicted  using  the  100,000-year  simulation  of  hurricanes.  To 
explore  this  comparison,  the  team  used  a  coarse  ADCIRC  grid  and  a  simple 
representation  of  the  hurricane  wind  field.  The  JPM  tracks  shown  in 
Figure  2-22  are  equally  spaced,  but  this  was  not  the  case  with  the  first  set  of 
156  hurricanes  (presented  in  Appendix  D).  Those  original  tracks  were 
obtained  by  removing  312  tracks  from  the  equally  spaced  set  of  468  hurri¬ 
canes.  This  first  set  of  156  storms  was  produced  during  a  project  review 
meeting  where  2/3  of  the  simulated  storms  were  removed  (those  corre¬ 
sponding  to  the  upper  and  lower  bound  B  values),  with  very  encouraging 
results  (i.e.,  small  changes  in  100-year  return  period  surge  elevations).  The 
removal  of  2/3  of  the  storms  results  in  unequally  spaced  storm  tracks,  and 
these  were  later  replaced  by  nominally  equally  spaced  tracks,  the  results  of 
which  are  presented  in  this  section  of  the  report.  The  results  of  the 
unequally  spaced  results  are  given  in  Appendix  D,  and  a  summary  of  the 
effect  of  storm  spacing  is  presented  at  the  conclusion  of  this  section. 

Figure  2-23  shows  the  locations  of  selected  sites  used  to  develop  plots  of  still 
water  level  vs.  return  period.  These  plots  are  used  to  compare  the  full  surge 
hazard  curves  arising  from  both  the  stochastic  and  JPM  simulations.  To 
develop  the  flood  hazard  curves  associated  with  the  JPM  hurricanes,  water 
elevations  were  computed  using  the  coarse  ADCIRC  grid  for  each  of  the 
156  JPM  simulated  storms.  Each  JPM  hurricane  has  a  probability  of 
occurrence  (conditional  on  a  hurricane  affecting  the  study  region)  equal  to 
the  product  of  the  individual  parameter  weights  given  in  Tables  2-1  through 
2-3.  The  sum  of  the  weights  over  all  three  storm  types  (North  Carolina 
landfalling,  Virginia/Delaware/New  Jersey  landfalling,  and  by-passing)  is 
equal  to  unity.  The  weights  in  Tables  2-1  through  2-3  have  to  be  multiplied 
by  0.202,  0.460  and  0.337,  respectively,  to  take  into  account  the  relative 
coastal/by-passing  segment  crossing  rates.  The  weights  are  further  multi¬ 
plied  by  1/3  to  take  into  account  that  three  values  of  RMW  are  used  for  each 
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Figure  2-22.  (concluded) 


pressure-heading-translation  speed  combination  (Vickery  and  Wadhera 
2008).  These  rates  are  further  multiplied  by  an  additional  term  equal  to 
i/Nb,  where  Nb  is  the  number  of  B  values  used  in  the  simulation  (Nb  is  one 
out  of  the  156  storm  cases  and  three  out  of  the  468  storm  cases). 

The  probability  that  a  hurricane-induced  water  elevation  is  exceeded 
during  time  period  t  is: 


Pt  (n  >  q0 )  =  1  -  (n  >  ho  I  x)  pt  (x)  (2-1) 

x=0 


where  P(ri>rio  \  x)  is  the  probability  that  the  water  elevation  77  is  greater 
than  t/o,  given  that  x  storms  occur,  and  pt(x)  is  the  probability  of  x  storms 
occurring  during  time  period  t.  From  Equation  2-1,  with  pt(x)  defined  as 
Poisson  and  defining  t  as  one  year,  the  annual  probability  of  exceeding  a 
storm  surge  elevation  is: 
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Figure  2-23.  Map  showing  the  locations  of  the  35  points  used  for  comparisons  of  storm  surge 
elevations  that  were  computed  using  the  full  storm  set  and  the  reduced  JPM  storm  set. 
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Pa(n>n0)  =  l-exp[-XP(Ti>Tlo)]  (2-2) 

where  A  represents  the  average  annual  number  of  storms  that  cross  any  of 
the  modeled  coastline  segments  and  P(Tj>rio)  is  the  probability  that  the 
water  elevation  rj  is  greater  than  rj(),  given  the  occurrence  of  any  one  storm. 
Again,  in  the  development  of  the  conditional  cumulative  distribution  func¬ 
tion  for  water  elevation,  P(r/>rio  \  x),  each  simulated  hurricane  used  to 
develop  the  distribution  has  a  probability  of  occurrence  of  Wi,  where  Wi  is 
the  product  of  the  individual  storm  weights  given  in  Tables  2-1  through  2-3 
(multiplied  by  the  appropriate  coastal,  RMW  and  B  weighting  factors).  The 
annual  occurrence  rate,  A,  is  0.1561. 

Figure  2-24  shows  the  comparisons  of  the  estimated  storm  surge  obtained 
using  the  JPM  set  of  storm  tracks  with  the  full  stochastic  set  for  a  return 
period  range  of  10  through  1,000  years.  The  storm  surge  estimates 
obtained  using  the  JPM  approach  are  reasonable  for  most  of  the  locations 
up  to  a  return  period  of  1,000  years.  The  cases  where  the  JPM  results  are 
notably  higher  than  the  stochastic  results  are  typically  at  the  end  of  bays 
and  estuaries.  Maps  showing  the  differences  between  the  two  sets  of  model 
results  on  a  node-by-node  basis  are  presented  later  in  this  report.  These 
maps  show  the  variation  in  the  differences  along  the  coast  of  the  Region 
III  study  area. 

Figure  2-25  presents  the  locations  (not  numbered)  for  all  near  coast  surge 
points  used  in  the  coarse  grid  analysis.  The  full  distribution  of  both  the 
relative  and  absolute  differences  between  the  stochastic  storm  set  and  the 
JPM  storm  set  are  presented  in  Figure  2-26  for  return  periods  of  10, 100 
and  500  years.  Figures  2-27  through  2-29  show  the  absolute  differences  in 
a  map  format.  Figures  2-30  through  2-32  present  the  same  information  as 
provided  in  Figures  2-27  through  2-29,  but  here  the  relative  differences 
are  presented  rather  than  the  absolute  differences.  Each  figure  presents 
the  maximum,  minimum,  and  mean  difference. 

Appendix  C  presents  the  results  for  468  equally  spaced  tracks  and  Appendix 
D  presents  the  results  for  156  equally  spaced  tracks.  Appendices  C  and  D 
also  provide  all  of  the  information  presented  in  Figures  2-20  through  2-32. 

Tables  2-4  and  2-5  summarize  the  maximum,  minimum,  mean,  RMS,  and 
the  5th  and  95th  percentile  differences  associated  with  the  data  presented  in 
Figures  2-24  through  2-32.  The  tables  also  present  the  difference  statistics 
derived  from  the  simulations  summarized  in  Appendices  C  and  D  for  the 


ERDC/CHL  TR-11-1;  Report  3 


33 


Figure  2-25.  Locations  of  all  storm  surge  points  used  to  derive  distributions  of  the  differences 
between  the  stochastic  storm  set  elevations  and  the  JPM  storm  set  elevations.  Storm  tracks 

are  shown  in  Figure  2-22. 

other  two  storm  sets  (Appendix  C  contains  the  468  storm  track  results  and 
Appendix  D  contains  the  156  equally  spaced  storms  track  results).  Table  2-4 
presents  the  statistics  for  the  differences  expressed  in  meters  and  Table  2-5 
presents  the  statistics  in  terms  of  the  percentage  differences.  In  all  cases,  a 
negative  difference  indicates  that  the  JPM  results  are  lower  than  those 
obtained  from  the  full  stochastic  simulation. 

The  results  indicate  that  the  468  storm  set  yields  the  least  amount  of 
variability  (as  defined  by  the  RMS  difference),  with  the  unequally  spaced 
reduced  storm  set  being  next  best.  In  the  case  of  the  100-year  return  period 
surge  elevations,  the  difference  range  as  defined  by  the  5th  to  95th  percentiles 
for  all  three  JPM  sets  is  less  0.26  m  (10.5  in.).  In  all  three  cases,  there  is  a 
negative  bias  in  the  estimates  of  the  10-year  return  period  surge  elevations, 
with  the  negative  bias  being  the  least  in  the  case  of  the  156  equally  spaced 
storm  tracks. 
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Figure  2-26.  Cumulative  distribution  functions  showing  the  differences  in  the  modeled  storm 
surge  elevations  at  all  nodal  points  shown  in  Figure  2-20.  The  differences  are  defined  as  the 
JPM  results  minus  the  stochastic  results.  Plots  on  the  left  present  the  relative  (percentage) 
difference  and  plots  on  the  right  show  the  absolute  difference,  expressed  as  meters. 
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Figure  2-27.  Difference  between  the  stochastic  and  JPM-predicted,  100-year  return  period 
storm  surge  elevations.  Positive  values  indicate  that  the  JPM  methodology  produces  higher 
storm  surge  values  than  the  stochastic  model.  Note  that  the  largest  positive  differences  occur 
on  the  north  side  of  Delaware  Bay,  along  the  New  Jersey  Coast. 

The  difference  statistics  for  the  500-year  return  period  JPM  storms  show  a 
small  high  bias,  with  the  largest  bias  evident  in  the  156  equally  spaced  storm 
set.  The  distribution  of  the  differences  has  a  long  positive  (overestimate)  tail 
that  is  not  evident  in  the  10  year  and  100  year  return  period  comparisons. 

Given  the  uncertainty  in  the  entire  modeling  process,  it  is  suggested  that 
the  equally  spaced  156  storms  presented  herein  be  used  as  the  simulation 
methodology  basis  for  developing  the  FIRMs. 
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Figure  2-28.  Difference  between  the  stochastic  and  JPM-predicted,  500-year  return  period 
storm  surge  elevations.  Positive  values  indicate  that  the  JPM  methodology  produces  higher 
storm  surge  values  than  the  stochastic  model.  Note  that  the  largest  positive  differences  occur 
on  the  north  side  of  Delaware  Bay,  and  at  the  most  extreme  inland  of  Chesapeake  Bay. 


Table  2-4.  Summary  statistics  of  differences  between  the  stochastic  and  JPM  storms  for  the 
three  different  JPM  storm  sets.  Differences  are  expressed  in  meters. 


156  Storms  -Equally 
Spaced  (Chapter  2) 

468  Storms  -Equally 
Spaced  (Appendix  C) 

156  Storms  -Unequally 
Spaced  (Appendix  D) 

10  yr 

100  yr 

500  yr 

10  yr 

100  yr 

500  yr 

10  yr 

100  yr 

500  yr 

Min 

-0.49 

-0.45 

-0.53 

-0.45 

-0.28 

-0.38 

-0.59 

-0.44 

-0.54 

Max 

0.04 

0.74 

1.00 

0.00 

0.52 

0.75 

0.01 

0.65 

0.81 

Mean 

-0.14 

-0.01 

0.18 

-0.19 

0.01 

0.20 

-0.26 

-0.01 

0.07 

RMS 

0.12 

0.16 

0.37 

0.08 

0.10 

0.26 

0.12 

0.14 

0.26 

5th 

-0.37 

-0.25 

-0.31 

-0.37 

-0.13 

-0.19 

-0.45 

-0.25 

-0.34 

95th 

0.00 

0.22 

0.82 

-0.07 

0.20 

0.65 

-0.07 

0.21 

0.53 
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Figure  2-29.  Difference  between  the  stochastic  and  JPM-predicted,  10-year  return  period 
storm  surge  elevations.  Positive  values  indicate  that  the  JPM  methodology  produces  higher 
storm  surge  values  than  the  stochastic  model.  There  is  a  clear  bias,  with  the  JPM  results 
being  lower  than  the  stochastic  results  for  all  locations. 


Table  2-5.  Summary  statistics  of  differences  between  the  stochastic  and  JPM  storms  for  the 
three  different  JPM  storm  sets.  Differences  are  expressed  as  a  percentage  of  the  stochastic 
storm  set  results  and  have  been  computed  on  a  point-by-point  basis. 


156  Storms  -Equally 
Spaced  (Chapter  2) 

468  Storms  -Equally 
Spaced  (Appendix  C) 

156  Storms  -Unequally 
Spaced  (Appendix  D) 

10  yr 

100  yr 

500  yr 

10  yr 

100  yr 

500  yr 

10  yr 

100  yr 

500  yr 

Min 

-73% 

-14% 

-16% 

-58% 

-12% 

-10% 

-68% 

-15% 

-14% 

Max 

11% 

45% 

22% 

0% 

26% 

20% 

3% 

32% 

17% 

Mean 

-27% 

0% 

4% 

-34% 

0% 

6% 

-46% 

0% 

2% 

RMS 

21% 

5% 

10% 

8% 

5% 

7% 

15% 

6% 

7% 

5th 

-60% 

-9% 

-11% 

-50% 

-9% 

-7% 

-63% 

-11% 

-10% 

95th 

0% 

7% 

19% 

-24% 

7% 

16% 

-18% 

8% 

12% 
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Figure  2-30.  Relative  (percentage)  difference  between  the  stochastic  and  JPM-predicted, 
100-year  return  period  storm  surge  elevations.  Positive  values  indicate  that  the  JPM  method¬ 
ology  produces  higher  storm  surge  values  than  the  stochastic  model.  Note  that  the  largest 
positive  differences  occur  on  the  north  side  of  Delaware  Bay,  along  the  New  Jersey  Coast. 
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Figure  2-31.  Relative  (percentage)  difference  between  the  stochastic  and  JPM- 
predicted,  500-year  return  period  storm  surge  elevations.  Positive  values 
indicate  that  the  JPM  methodology  produces  higher  storm  surge  values  than  the 
stochastic  model.  Note  that  the  largest  positive  differences  occur  on  the  north 
side  of  Delaware  Bay,  along  the  New  Jersey  Coast,  and  at  the  inland  most 
extreme  of  Chesapeake  Bay. 
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Figure  2-32.  Relative  (percentage)  difference  between  the  stochastic  and  JPM-predicted, 
10-year  return  period  storm  surge  elevations.  Positive  values  indicate  that  the  JPM  methodology 
produces  higher  storm  surge  values  than  the  stochastic  model.  There  is  a  clear  bias  with  the 
JPM  results  being  lower  than  the  stochastic  results  for  all  locations. 
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3  Extratropical  Cyclones 

The  selection  and  development  of  extratropical  cyclones  (ETC)  and 
associated  wind  and  pressure  fields  for  the  FEMA  Region  III  storm  surge 
study  was  performed  by  Oceanweather,  Inc.  Initially,  a  set  of  30  extratropical 
cyclones  were  identified  for  hindcast  based  on  a  ranking  of  historical 
measured  water  levels  from  seven  stations  in  the  NOS  archive.  After  the 
completion  of  this  work,  in  November  2009,  ETC  Ida  made  a  significant 
impact  on  the  study  area.  Ida  originated  as  Hurricane  Ida  in  the  Gulf  of 
Mexico  and  resulted  in  water  levels  8. 0-9.0  ft  above  MLLW  at  some  locations 
in  the  study  domain  and  5.0-6.0  ft  above  MLLW  in  the  Chesapeake  Bay. 
Sample  ETC  Ida  water  levels  from  two  NOS  stations  in  the  Chesapeake  Bay 
appear  in  Figure  3-1.  Predicted  (tides  only),  observed,  and  residual  water 
levels  are  depicted.  The  significant  increase  of  observed  water  levels  above  the 
tidal  prediction  is  a  measure  of  the  Ida  storm  surge  and  flooding  at  this 
location.  As  a  result  of  this  widespread  impact,  Ida  was  added  to  the  list  of 
extratropical  cylones  selected  for  the  study.  Furthermore,  Ida  has  been  added 
to  the  list  of  validation  storms  that  will  be  used  to  quantify  the  performance  of 
the  modeling  system  prior  to  production. 

The  documentation  for  this  work  was  provided  by  Oceanweather  and 
appears  in  a  series  of  two  appendices.  Appendix  A,  updated  after  ETC  Ida 
was  added,  provides  a  summary  of  the  work  performed  by  Oceanweather 
to  identify  the  extratropical  cyclone  dataset  and  reconstruct  these  events. 
Appendix  B  describes  the  extratropical  cyclone  selection  process  for  the 
study.  Included  are  the  NOS  station  data  used  to  rank  the  events,  and 
samples  of  the  reconstructed  water  levels  and  pressure  fields. 
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Figure  3-1.  NOS  water  levels  at  Chesapeake  Bay  Bridge  Tunnel  (upper)  and  Sewells  Point,  VA 
(lower)  during  Extratropical  Cyclone  Ida  in  November  2009. 
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FEMA  Region  III  Storm  Surge  Modeling 
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Cos  Cob,  CT,  USA 
Tel:  203-661-3091 
Fax:  203-661-6809 
Email:  oceanwx@oceanweather.com 
Web:  www.oceanweather.com 

Introduction 

This  report  describes  work  performed  under  contract  W912BU-08-P-0291 
for  the  US  Army  Engineer  District,  Philadelphia,  for  storm  reanalysis  in 
support  of  FEMA  Region  III  storm  surge  modeling.  The  objective  of  the 
project  is  to  update  the  coastal  flood  risk  analysis  for  Chesapeake  Bay  (CB) 
and  eventually  all  of  Region  III.  The  project  will  incorporate  the  latest 
approaches,  modeling,  and  statistical  methodologies  applied  in  similar 
studies  carried  out  post- Katrina  to  Regions  VI  (Gulf  of  Mexico)  and  IV 
(Southeast  Coast  of  US).  Region  III  is  subject  to  both  extreme  extratropical 
cyclones  (ETC)  and  tropical  cyclones  (TC),  so  both  regimes  need  to  be 
addressed. 

This  project  consists  of  the  following  major  tasks  and  deliverables: 

1.  ETC  storm  selection 

2.  Hindcast  of  25  ETC  storms  for  Chesapeake  Bay  support 

3.  Hindcast  of  5.0  additional  ETC  and  2.0  TC  for  Region  III  support 
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4.  Transfer  of  TC96  Tropical  Planetary  Boundary  Layer  model 

5.  Additional  ETC  case  of  the  “Nor’Ida”  storm  of  November  2009 

Extratropical  storm  selection 

A  total  of  30  ETCs  were  identified  for  hindcast  using  historical  measured 
water  levels  from  the  National  Ocean  Service  (NOS)  archive.  Seven  stations 
in  the  Region  III  area  with  long  term  continuous  records  during  the  1975- 
2008  period  were  identified  for  the  storm  selection  process.  A  summary 
report  titled  Extra-Tropical  Storm  Selection  In  Support  ofFEMA  Region 
III  Storm  Surge  Modeling  provides  details  about  the  selection  process. 

After  the  30  ETC  storm  work  was  completed,  an  extratropical  system  (so- 
called  “Nor’Ida”)  which  originated  as  Hurricane  Ida  in  the  Gulf  of  Mexico 
impacted  the  study  area  during  the  period  of  November  11-15,  2009.  This 
storm  resulted  in  water  levels  8. 0-9.0  ft  above  MLLW  at  The  Battery,  New 
York  (NOS  8518750),  Sandy  Hook,  New  Jersey  (NOS  8531680),  and  Cape 
May,  New  Jersey  (8536110),  and  impacted  numerous  other  NOS  gauges  in 
the  study  area.  This  storm  was  hindcast  using  the  same  ETC  methodology 
outlined  below  and  added  to  the  project. 

Extratropical  storm  hindcast 

The  30  ETC  storms  plus  the  Nor’Ida  storm  identified  for  hindcast  are  shown 
in  Table  Ai  and  include  the  approximately  five-day  spin-up  and  three-day 
spin-down  period  from  the  storm’s  peak  surge.  Storms  were  hindcast  using 
the  Interactive  Objective  Kinematic  Analysis  (IOKA)  system  as  described  in 
Cox  et  al.  1995.  Wind  and  pressure  data  were  obtained  from  the  following 
sources: 

1.  National  Centers  for  Environmental  Prediction/National  Center  for 
Atmospheric  Research  (NCEP/NCAR)  reanalysis  wind  and  pressure  fields 

2.  Buoy  and  Coastal  Manned  Station  (C-MAN)  data  from  National  Data  Buoy 
Center  (NDBC) 

3.  National  Weather  Service  (NWS)  and  National  Ocean  Service  (NOS)  land 
station  data  provided  by  the  National  Climatic  Data  Center  (NCDC) 

4.  Ship  reports  from  Comprehensive  Ocean- Atmosphere  Data  Set  (COADS) 

5.  Scatterometer  wind  estimates  from  the  QUIKSCAT,  ERS-i,  ERS-2 
instruments 
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Table  Al.  Extratropical  storm  hindcast  list. 


Storm  Reference 

Start  (UTC) 

End  (UTC) 

19750926 

Sep-21-1975  00:00 

Sep-29-1975  00:00 

19780126 

Jan-21-1978  00:00 

Jan-29-1978  00:00 

19780427 

Apr-22-1978  00:00 

Apr-30-1978  00:00 

19790226 

Feb-21-1979  00:00 

Mar-01-1979  00:00 

19821025 

0ct-20-1982  00:00 

Oct-28-1982  00:00 

19830404 

Mar-30-1983  00:00 

Apr-07-1983  00:00 

19840216 

Feb-11-1984  00:00 

Feb-19-1984  00:00 

19840329 

Mar-24-1984  00:00 

Apr-01-1984  00:00 

19850212 

Feb-07-1985  00:00 

Feb- 15-1985  00:00 

19851107 

Nov-02-1985  00:00 

Nov- 10-1985  00:00 

19870102 

Dec-28- 1986  00:00 

Jan-05-1987  00:00 

19911031 

Oct-26-1991 00:00 

Nov-03-1991  00:00 

19920105 

Dec-31-1991  00:00 

Jan-08- 1992  00:00 

19921211 

Dec-06-1992  00:00 

Dec-14-1992  00:00 

19930305 

Feb-28-1993  00:00 

Mar-08-1993  00:00 

19931128 

Nov-23-1993  00:00 

Dec-01-1993  00:00 

19940104 

Dec-30-1993  00:00 

Jan-07-1994  00:00 

19960108 

Jan-03-1996  00:00 

Jan-11-1996  00:00 

19960119 

Jan-14-1996  00:00 

Jan-22-1996  00:00 

19961009 

0ct-04-1996  00:00 

Oct-12-1996  00:00 

19980128 

Jan-23-1998  00:00 

Jan-31-1998  00:00 

19980205 

Jan-31-1998  00:00 

Feb-08-1998  00:00 

19980309 

Mar-04-1998  00:00 

Mar-12-1998  00:00 

20000125 

Jan-20-2000  00:00 

Jan-28-2000  00:00 

20021017 

0ct-12-2002  00:00 

0ct-20-2002  00:00 

20051025 

0ct-20-2005  00:00 

0ct-28-2005  00:00 

20061007 

0ct-02-2006  00:00 

Oct- 10-2006  00:00 

20061028 

0ct-23-2006  00:00 

0ct-31-2006  00:00 

20061122 

Nov- 17-2006  00:00 

Nov-25-2006  00:00 

20080512 

May-07-2008  00:00 

May-15-2008  00:00 

20091113 

Nov-08-2009  00:00 

Nov- 16-2 009  00:00 
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Wind  fields  from  individual  storms  were  reanalyzed  using  kinematic 
analysis,  a  man-intensive  process  for  which  wind  speeds  (isotachs)  and 
wind  directions  (streamlines)  are  hand-drawn  to  best  represent  the 
available  observations  while  preserving  the  primary  meteorological 
principles  of  storm  development  and  continuity.  Figure  Ai  shows  an 
example  of  a  kinematic  analysis  valid  during  the  October  2006  event. 
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Figure  Al.  Kinematic  analysis  of  wind  speed  (red  knots)  and  wind  direction  (streamlines, 
black)  valid  October  7,  2006  00:00  UTC. 
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Wind  fields  were  developed  on  two  grid  systems:  a  basin  grid  covered  the 
domain  15-58N,  85-  55 W  at  a  0.125-deg  latitude-longitude  grid  with  time 
step  of  30  min.  A  fine  grid  covered  the  domain  36-40N,  78-74W  at 
0.025-deg  with  15  min  time  step.  All  winds  are  10-m  neutral  representing 
a  30-min  average.  Pressures  from  the  NCEP/NCAR  data  were  interpolated 
onto  the  same  working  grids  as  the  wind  fields. 

Tropical  cyclone  hindcasts 

The  reference  set  of  historical  storms  selected  for  calibration  of  the 
modeling  consist  of  Isabel  (2003)  and  Ernesto  (2006).  Copious 
meteorological  data  exist  and  have  been  processed  to  develop  wind  and 
pressure  fields  to  drive  ocean  response  models.  The  pertinent  data  sets  for 
wind  field  analysis  consist  of 

1.  historical  track  and  intensity  data  from  NOAA’s  (National  Oceanographic 
and  Atmospheric  Administration)  HURD  AT  (Hurricane  Data)  archive; 

2.  aircraft  reconnaissance  obtained  from  NOAA  and  US  Air  Force  hurricane 
hunter  aircraft,  including  vortex  messages  and  continuous  flight  level  wind 
speed,  direction,  DValue; 

3.  gridded  and  image  fields  of  marine  surface  wind  composites  from  the 
HRD  (Hurricane  Research  Division)  HWnd  analysis; 

4.  synoptic  observations  from  NOAA  buoy  and  C-MAN  (Coastal  Manned) 
stations; 

5.  synoptic  observations  from  transient  ships,  coastal  and  land  stations; 

6.  composite  NWS  (National  Weather  Service)  radar  imagery; 

7.  loops  of  NOAA  visual,  infrared  and  water  vapor  imagery; 

8.  NCEP  (National  Center  for  Environmental  Prediction)  model  wind  fields; 

9.  QUIKSCAT  scatterometer  winds; 

10.  TOPEX  altimeter  winds  and  waves; 

11.  ERS-2  altimeter  winds  and  waves;  and 

12.  data  from  offshore  platforms  equipped  with  met  packages. 

The  meteorological  parameters  of  greatest  importance  are  the  hurricane's 
track,  forward  speed,  radius  of  maximum  wind,  barometric  pressure 
distribution  from  the  hurricane  center  to  the  periphery,  and  the  ambient 
pressure  system  through  which  the  hurricane  is  moving.  In  the  hindcast 
method,  these  parameters  are  used  to  drive  a  numerical  model  (Thompson 
and  Cardone  1996)  (TC96)  to  determine  a  first  guess  of  the  surface  wind 
field  using  OWI,  Inc.’s  (OWI’s)  tropical  analyst  workstation  (Cox  and 
Cardone  2000);  next,  the  available  direct  wind  measurements  were 


ERDC/CHL  TR-11-1;  Report  3 


50 


assimilated  into  the  first  guess  fields  to  provide  a  refined  specification  using 
OWI’s  Interactive  Objective  Kinematic  Analysis  (IOKA)  system  (Cox  et  al. 
1995).  Since  the  1950s,  daily  reconnaissance  flights  have  been  made  into 
most  storms  by  Air  Force  and  NOAA  personnel  to  collect  pressure,  tempera¬ 
ture,  and  upper  level  wind  data.  Eye  shape  and  diameter  observations  are 
also  recorded.  These  data  are  used  to  fit  the  snapshot  parameters  required 
to  drive  the  TC96  model. 

Wind  and  pressure  fields  were  developed  on  the  same  basin/fine  grid 
systems  applied  in  the  ETC  hindcast  work.  Pressure  output  from  the  TC96 
model  within  the  core  of  the  system  was  blended  into  the  larger  domain  of 
the  NCEP/NCAR  reanalysis  pressures. 

Validation  of  Ernesto  (2006_06)  and  Extratropical  cyclone  20051025 

Table  A2  and  Figure  A2  show  the  closest  buoys  and  C-MAN  station  data 
locations  for  which  data  were  verified  against  the  hindcast.  All  data  were 
obtained  from  the  National  Data  Buoy  Center  and  adjusted  for  both  height 
and  stability  to  a  10-m  reference  level.  Where  available,  three  consecutive 
10-min  wind  observations  were  averaged  1-1-1  to  produce  a  30-min  mean 
wind  for  direct  comparison  to  the  hindcast  sets. 

Time  histories  of  each  station  and  storm  are  shown  in  Figures  A3 
(20051025)  and  A4  (Ernesto).  Mean  wind  speed  bias  during  20051025 
was  0.22  m/sec  with  correlation  coefficient  of  92  %,  wind  direction  bias 
was  under  1.0  deg.  Mean  wind  bias  during  Ernesto  was  0.32  m/sec  with 
correlation  coefficient  of  98  %;  wind  direction  bias  was  also  within  1.0  deg. 


Table  A-2  Buoys  and  C-MAN  validation  stations. 


Station 

Latitude 

Longitude 

44009 

38.464 

-74.702 

44014 

36.611 

-74.836 

CHLV2 

39.910 

-75.710 

TPLM2 

38.898 

-76.437 
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Figure  A2.  Verification  locations  (image  courtesy  of  NDBC). 


Figure  A3.  Comparison  of  wind  speed  (top)  and  direction  (bottom)  for  four  locations  during  extra- 
tropical  cyclone  20051025.  Red  are  NDBC  observations,  black  is  the  hindcast. 


ERDC/CHL  TR-11-1;  Report  3 


52 


Figure  A4.  Comparison  of  wind  speed  (top)  and  direction  (bottom)  for  four  locations  during  Ernesto 
(2006_06).  Red  are  NDBC  observations,  black  is  the  hindcast. 


TC96  tropical  planetary  boundary  layer  model  transfer 

The  source  code  and  documentation  for  the  TC96  tropical  boundary  layer 
model  use  application  were  provided  to  RENCI  for  application  in  the 
Region  III  modeling  work.  The  model  is  contained  within  a  set  of 
FORTRAN  modules,  which  was  delivered  along  with  documentation 
pertaining  to  the  model  implementation. 
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Approach 

This  report  describes  the  extratropical  cyclone  selection  process  for  the 
FEMA  Region  III  coastal  storm  surge  study.  The  primary  area  addressed 
includes  both  Chesapeake  Bay  and  Delaware  Bay  and  applies  storm-ranking 
of  surge  events  based  on  NOS  measurements.  Applying  the  measured  NOS 
water  level  data  —  assuming  a  long  continuous  record  is  available  at 
multiple  stations  —  is  consistent  with  other  recent  FEMA  studies,  such  as 
that  which  was  performed  in  North  Carolina.  A  total  of  30  storm  events 
were  identified  for  hindcast. 
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NOS  station  data 

Historical  NOS  water  level  data  are  available  for  download  from  the  NOS 
website  and  include  both  a  measured  and  predicted  (tide)  component.  The 
residual  (measured-tide)  water  level  was  computed  for  all  stations  to 
remove  the  tidal  response  so  that  the  storm  selections  are  based  on  forcing, 
rather  than  tidal  considerations.  Measured  data  in  the  Chesapeake  and 
Delaware  Bays  are  shown  in  Figure  Bi  and  include  20  measurement 
locations  with  indicated  periods  of  record  of  at  least  15-years  length.  In 
order  to  keep  the  storm  selection  process  as  unbiased  as  possible,  multiple 
long-term  records  of  water  level  data  are  required.  Stations  with  large  gaps, 
or  short  time  records  were  not  considered.  Figure  B2  shows  the  measured 
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Figure  Bl.  Available  NOS  station  data  in  Chesapeake  and  Delaware  Bay  with  records  greater 

than  15  years  length. 
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Figure  B2.  Historical  residual  water  level  measurements  at  Chesapeake  City,  MD  (NOS 
8573927)  shown  as  an  example  of  a  record  with  large  gaps. 


residual  water  level  at  Chesapeake  City,  Maryland  (NOS  station  8573927), 
for  which  there  are  no  reports  from  1985-2004.  When  each  of  the  NOS 
records  are  examined,  only  seven  stations  (Figure  B3  and  Table  Bi)  are 
found  to  have  good  coverage  during  the  period  January  26, 1975  (based  on 
NOS  8638863)  to  August  31,  2008  (present  time  when  the  selection  was 
performed). 

Storm  ranking 

At  each  selected  NOS  station,  the  water  level  peaks  greater  than  the  99th 
percentile  residual  water  level  were  identified  individually  as  candidate 
storm  events.  Figure  B4  shows  the  residual  water  level  record  at  Chesapeake 
Bay  Bridge  Tunnel,  VA  (NOS  8638863)  where  peaks  are  identified 
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Figure  B3.  Selected  NOS  stations  applied  in  the  storm  selection  process. 


Table  Bl.  Select  NOS  stations  applied  in  the  storm  selection  process. 


Station  ID 

Lat 

Long 

Record 

Began 

(YMD) 

Record 

Ended 

(YMD) 

Station  Location 

8638863 

36.966 

-76.113 

19750126 

20080831 

Chesapeake  Bay  Bridge,  VA 

8638610 

36.946 

-76.33 

19270701 

20080831 

Sewells  Point,  VA 

8575512 

38.983 

-76.466 

19280801 

20080731 

Annapolis,  MD 

8574680 

39.266 

-76.578 

19020701 

20080831 

Baltimore,  MD 

8557380 

38.781 

-75.12 

19570101 

20080831 

Lewes,  DE 

8536110 

38.968 

-74.96 

19651101 

20080831 

Cape  May,  NJ 

8594900 

38.873 

-77.021 

19310401 

20080831 

Washington,  DC 
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Figure  B4.  Peaks  (red)  greater  that  the  99th  percentile  residual  water  level  identified  at 
Chesapeake  Bay  Bridge  Tunnel,  VA  NOS  station. 


(filled  red  circles  in  Figure  B4)  in  the  January  1975  to  August  2008  period. 
Storm  peaks  were  then  ranked  by  station  and  the  top  ranked  events 
matched  with  those  from  the  other  NOS  stations  within  a  3.0-day  window. 

A  population  of  30  events  was  obtained  by  taking  the  top-ranked  storms 
from  each  of  the  NOS  stations  sequentially  and  removing  dates  found  to  be 
tropical  in  nature,  as  determined  by  historical  tropical  storm/hurricane 
tracks.  During  this  period,  the  tropical  cyclones  of  David  (1979),  Gloria 
(1985),  Josephine  (1996),  Isabel  (2003),  and  Ernesto  (2006)  were  identified 
and  removed  from  the  storm  list.  Multiple  storm  peaks  at  a  single  NOS  site 
within  a  3.0-day  period  were  merged  into  a  single  storm  event.  The  storm- 
surge-ranked  table  represents  at  least  the  top  six  measured  events  at  each 
individual  station  (Table  B2).  Storms  with  missing  surge  values  typically 
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indicate  that  the  measured  water  level  was  under  the  99th  percentile 
threshold,  rather  than  missing  from  the  archive  entirely. 

NOS  stations  in  the  mid-Chesapeake  region  and  between  the  Chesapeake 
and  Delaware  Bays  were  removed  in  the  selection  process  due  to  poor 
temporal  coverage.  However,  their  removal  leaves  a  large  portion  of  the 
mid-Chesapeake  not  represented  by  a  NOS  measurement  station  for  the 
storm  ranking.  To  test  for  missing  storms,  a  separate  storm  ranking  was 
performed  for  the  period  1998-2005  using  the  Lewisetta,  Virginia  (NOS 
8635750),  Wachapreauge,  Virginia  (NOS  8631044);  and  Brandywine  Shoal 
Light,  Delaware  (NOS  8555889)  stations.  Table  B3  shows  the  top  five  extra- 
tropical  events  at  each  station.  This  ranking  indicates  that  four  of  the  five 
top  storms  at  Lewistetta,  three  of  five  at  Wachapreauge,  and  three  of  five  at 
Brandywine  Shoal  Light  were  identified  for  hindcast,  and  the  top  events  in 
the  1998-2003  period  at  each  station  are  also  represented  in  the  hindcast. 
Figure  B5  shows  plots  of  residual  water  level  and  sea  level  pressures  for 
each  storm. 
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Table  B2.  Extratropical  cyclones  selected  from  the  period  Jan-1975  to  Aug-2008  based  on  seven  NOS  stations. 
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Table  B3. 1998-2005  storms  found  at  mid-Chesapeake  NOS  stations. 
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Figure  B5.  (Sheet  2  of  60) 
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Figure  B5.  (Sheet  3  of  60) 
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Figure  B5.  (Sheet  6  of  60) 


ERDC/CHL  TR-11-1;  Report  3 


67 
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Figure  B5.  (Sheet  8  of  60) 
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Figure  B5.  (Sheet  14  of  60) 


^1 

Ol 
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Figure  B5.  (Sheet  22  of  60) 


ERDC/CHL  TR-11-1;  Report  3 


83 
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Figure  B5.  (Sheet  50  of  60) 


Figure  B5.  (Sheet  51  of  60) 


ERDC/CHL  TR-11-1;  Report  3  111 


ERDC/CHL  TR-11-1;  Report  3 


112 
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Figure  B5.  (Sheet  57  of  60) 


ERDC/CHL  TR-11-1;  Report  3 


118 
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Appendix  C:  468  JPM  Track  Results 

This  Appendix  discusses  the  results  of  and  comparisons  to  the  stochastic 
model  values  resulting  from  the  468  storm  JPM  set. 

The  probability  distributions  of  the  key  landfall  parameters  and  associated 
weights  used  in  the  JPM  approach  using  the  468  storms  are  the  same  as 
those  used  for  the  156-storm  case  and  are  reflected  in  Figures  2-14  through 
2-17  and  Tables  2-1  through  2-3,  as  given  in  Chapter  2. 

Figures  Cl  and  C2  represent  comparisons  of  the  distributions  of  the 
landfall  values  of  central  pressure,  etc.  from  the  468  storm  JPM  set  to  the 
stochastic  simulation  results.  The  comparisons  are  essentially  the  same  as 
those  discussed  in  Chapter  2  for  the  156-storm  case,  with  the  exception 
that  the  complete  distribution  of  B  is  retained  in  this  larger  set  that 
contains  both  the  mean  model  for  B,  as  well  as  the  upper  and  lower  bound 
models. 

Figure  C3  presents  the  nine  sets  of  tracks  used  in  the  468-JPM  storm  set. 
Figure  C4  shows  the  comparisons  of  the  estimated  storm  surge  obtained 
using  the  JPM  set  of  storm  tracks  and  full  stochastic  set  for  a  return  period 
range  of  10  through  1,000  years.  As  seen  in  Chapter  2,  the  storm  surge 
estimates  obtained  using  the  JPM  approach  are  reasonable  for  most  of  the 
locations  up  to  a  return  period  of  1,000  years.  The  cases  where  the  JPM 
results  are  notably  higher  than  the  stochastic  results  are  typically  at  the 
end  of  bays  and  estuaries. 

The  full  distribution  of  both  the  relative  and  absolute  differences  between 
the  stochastic  and  the  JPM  storm  sets  are  presented  in  Figure  C5  for 
return  periods  of  10, 100  and  500  years.  Figures  C6  through  C8  represent 
the  absolute  differences  in  a  map  format.  Figures  C9  through  C11  present 
the  same  information  as  given  in  Figures  C6  through  C8,  but  here  the 
differences  are  presented  in  a  relative  sense  rather  than  the  absolute 
errors.  Each  figure  presents  the  maximum,  minimum,  and  mean  error. 
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Figure  Cl.  Comparison  of  JPM  and  stochastic  representations  of  hurricane  parameters  at 
landfall  along  the  VA/DE/NJ  coastal  segment.  Heading  is  defined  as  the  direction  of  motion  of 
the  storm  measured  clockwise  from  north.  Storms  with  a  positive  heading  have  an  eastward 
component.  The  translation  speed  is  denoted  c,  the  central  pressure  difference  is  denoted  dP 
and  RMW  is  the  radius  to  maximum  winds. 
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Figure  C2.  Comparison  of  JPM  and  stochastic  representations  of  hurricane  parameters  at 
landfall  along  the  NC  coastal  segment.  Heading  is  defined  as  the  direction  of  motion  of  the 
storm  measured  clockwise  from  north.  Storms  with  a  positive  heading  have  an  eastward 
component.  The  translation  speed  is  denoted  c,  the  central  pressure  difference  is  denoted  dP 
and  RMW  is  the  radius  to  maximum  winds. 
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Figure  C4.  Comparisons  of  storm  surge  elevations  computed  using  the  stochastic  and  the 
JPM  storm  sets  represented  with  468  simulated  tracks. 
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Figure  C5.  Cumulative  distribution  functions  showing  the  differences  in  the  modeled  storm 
surge  elevations  at  all  nodal  points  shown  in  Figure  2-20.  The  differences  are  defined  as  the 
JPM  results  minus  the  stochastic  results.  Plots  on  the  left  present  the  relative  difference  and 
plots  on  the  right  show  the  absolute  difference,  expressed  as  meters. 
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Figure  C6.  Difference  between  the  stochastic  and  JPM-predicted  100-year  return  period  storm 
surge  elevations.  Positive  values  indicate  that  the  JPM  methodology  produces  higher  storm 
surge  values  than  the  stochastic  model.  The  largest  positive  differences  occur  on  the  north 
side  of  Delaware  Bay,  along  the  New  Jersey  coast. 
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Figure  C7.  Difference  between  the  stochastic  and  JPM-predicted  500-year  return  period  storm 
surge  elevations.  Positive  values  indicate  that  the  JPM  methodology  produces  higher  storm 
surge  values  than  the  stochastic  model.  Note  the  largest  positive  differences  occur  on  the 
north  side  of  Delaware  Bay,  along  the  New  Jersey  coast  and  at  the  inland  most  extreme  of  the 

Delaware  Bay. 
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Figure  C8.  Difference  between  the  stochastic  and  JPM-predicted  10-year  return  period  storm 
surge  elevations.  Positive  values  indicate  that  the  JPM  methodology  produces  higher  storm 
surge  values  than  the  stochastic  model.  There  is  a  clear  bias,  with  the  JPM  results  being  lower 
than  the  stochastic  results  for  all  locations. 
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Figure  C9.  Relative  (percentage)  difference  between  the  stochastic  and  JPM-predicted 
100-year  return  period  storm  surge  elevations.  Positive  values  indicate  that  the  JPM 
methodology  produces  higher  storm  surge  values  than  the  stochastic  model.  The  largest 
positive  differences  occur  on  the  north  side  of  Delaware  Bay,  along  the  New  Jersey  coast. 
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Figure  CIO.  Relative  (percentage)  difference  between  the  stochastic  and  JPM-predicted 
500-year  return  period  storm  surge  elevations.  Positive  values  indicate  that  the  JPM 
methodology  produces  higher  storm  surge  values  than  the  stochastic  model.  The  largest 
positive  differences  occur  on  the  north  side  of  Delaware  Bay,  along  the  New  Jersey  coast  and 
at  the  inland  most  extreme  of  the  Delaware  Bay. 
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Figure  Cll.  Relative  (percentage)  difference  between  the  stochastic  and  JPM-predicted 
10-year  return  period  storm  surge  elevations.  Positive  values  indicate  that  the  JPM 
methodology  produces  higher  storm  surge  values  than  the  stochastic  model.  There  is  a  clear 
bias,  with  the  JPM  results  being  lower  than  the  stochastic  results  for  all  locations. 
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Appendix  D:  156  Unequally  Spaced  JPM  Track 
Results 

This  Appendix  discusses  the  results  and  comparisons  to  the  stochastic 
model  values  resulting  from  the  156  unequally  spaced  JPM  storm  set.  This 
set  of  tracks  was  developed  by  removing  2/3  of  the  tracks  from  the  full  set 
of  468,  which  included  the  modeling  of  the  variance  associated  with  the 
Holland  B  parameter.  The  removal  of  2/3  of  the  storms  resulted  in 
unequally  spaced  tracks,  but  it  enabled  the  rapid  assessment  of  the  impact 
of  removing  the  results  associated  with  the  upper  and  lower  bound 
representations  of  the  Holland  B  parameter  without  the  need  to  rerun  the 
ADCIRC  simulations. 

The  probability  distributions  of  the  key  landfall  parameters  and  associated 
weights  used  in  the  JPM  are  the  same  as  those  used  for  the  156-storm  case 
and  are  reflected  in  Figures  2-14  through  2-17  and  Tables  2-1  through  2-3, 
(Chapter  2)  Therefore  they  will  not  be  presented  here  again. 

Comparisons  of  the  distributions  of  the  landfall  values  of  central  pressure, 
etc.  between  the  JPM  set  to  the  stochastic  simulation  results  are  the  same 
as  those  seen  in  the  156  set  of  equally  spaced  tracks  presented  in  Section 
2.5  and  are  not  reproduced  here. 

Figure  Di  presents  the  nine  sets  of  tracks  used  in  this  storm  set  of  156 
equally  spaced  tracks.  Figure  D2  shows  the  comparisons  of  the  storm 
surge  estimated  using  the  JPM  set  of  storm  tracks  and  full  stochastic  set 
for  a  return  period  range  of  10  through  1,000  years.  As  seen  in  Chapter  2, 
Appendix  C,  and  again  here,  the  storm  surge  estimates  calculated  using 
the  JPM  approach  are  reasonable  for  most  of  the  locations,  up  to  a  return 
period  of  1,000  years.  The  cases  where  the  JPM  results  are  notably  higher 
than  the  stochastic  results  are  typically  at  the  end  of  bays  and  estuaries. 

The  full  distribution  of  both  the  relative  and  absolute  differences  between 
the  stochastic  and  the  JPM  storm  sets  are  presented  in  Figure  D3  for 
return  periods  of  10, 100,  and  500  years.  Figures  D4  through  D6  show  the 
absolute  differences  in  a  map  format.  Figures  D7  through  D9  present  the 
same  information  as  given  in  Figures  D4  through  D6,  but  here  the 
differences  are  presented  in  a  relative  sense  rather  than  as  absolute  errors. 
Each  figure  presents  the  maximum,  minimum,  and  mean  error. 
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Figure  D2  Comparisons  of  storm  surge  elevations  computed  using  the  stochastic  and  the  JPM 
storm  sets,  represented  with  156  unequally  spaced  simulated  tracks. 
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Figure  D3.  Cumulative  distribution  functions  showing  the  differences  in  the  modeled  storm 
surge  elevations  at  all  nodal  points  shown  in  Figure  2-20.  The  differences  are  defined  as  the 
JPM  results  minus  the  stochastic  results.  Plots  on  the  left  represent  the  relative  difference 
and  plots  on  the  right  show  the  absolute  difference,  expressed  as  meters. 
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Figure  D4.  Difference  between  the  stochastic-  and  JPM-predicted  100-year  return  period 
storm  surge  elevations.  Positive  values  indicate  that  the  JPM  methodology  produces  higher 
storm  surge  values  than  the  stochastic  model.  The  largest  positive  differences  occur  on  the 
north  side  of  Delaware  Bay,  along  the  New  Jersey  coast. 


ERDC/CHL  TR-11-1;  Report  3 


140 


Figure  D5.  Difference  between  the  stochastic-  and  JPM-predicted  500-year  return  period 
storm  surge  elevations.  Positive  values  indicate  that  the  JPM  methodology  produces  higher 
storm  surge  values  than  the  stochastic  model.  The  largest  positive  differences  occur  on  the 
north  side  of  Delaware  Bay,  along  the  New  Jersey  coast  and  at  the  inland  most  extreme  of  the 

Delaware  Bay. 
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Figure  D6.  Difference  between  the  stochastic-  and  JPM-predicted  10-year  return  period  storm 
surge  elevations.  Positive  values  indicate  that  the  JPM  methodology  produces  higher  storm 
surge  values  than  the  stochastic  model.  There  is  a  clear  bias,  with  the  JPM  results  being  lower 
than  the  stochastic  results  for  all  locations. 
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Figure  D7.  Relative  (percentage)  difference  between  the  stochastic-  and  JPM-predicted 
100-year  return  period  storm  surge  elevations.  Positive  values  indicate  that  the  JPM 
methodology  produces  higher  storm  surge  values  than  the  stochastic  model.  The  largest 
positive  differences  occur  on  the  north  side  of  Delaware  Bay,  along  the  New  Jersey  coast. 
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Figure  D8.  Relative  (percentage)  difference  between  the  stochastic-  and  JPM-predicted 
500-year  return  period  storm  surge  elevations.  Positive  values  indicate  that  the  JPM 
methodology  produces  higher  storm  surge  values  than  the  stochastic  model.  The  largest 
positive  differences  occur  on  the  north  side  of  Delaware  Bay,  along  the  New  Jersey  coast  and 
at  the  inland  most  extreme  of  the  Delaware  Bay. 
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Figure  D9.  Relative  (percentage)  difference  between  the  stochastic-  and  JPM-predicted 
10-year  return  period  storm  surge  elevations.  Positive  values  indicate  that  the  JPM 
methodology  produces  higher  storm  surge  values  than  the  stochastic  model.  There  is  a  clear 
bias  with  the  JPM  results  being  lower  than  the  stochastic  results  for  all  locations. 
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